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Abstract 



We demonstate the existence of a large N phase transition with respect to the 't 
Hooft coupling in q-deformed Yang-Mills theory on S 2 . The strong coupling phase is 
characterized by the formation of a clump of eigenvalues in the associated matrix model 
of the Douglas-Kazakov (DK) type pQ. By understanding this in terms of instanton 
contributions to the q-deformed Yang- Mills theory, we gain some insight into the strong 
coupling phase as well as probe the phase diagram at nonzero values of the 6 angle. The 
^ ■ Ooguri-Strominger-Vafa [2] relation of this theory to topological strings on the local 

Calabi-Yau 0{— p) ® 0{p — 2) — ► P 1 via a chiral decompostion at large N 3 a , motivates 
us to investigate the phase structure of the trivial chiral block, which corresponds to 
the topological string partition function, for p > 2. We find a phase transition at a 
different value of the coupling than in the full theory, indicating the likely presence of 
a rich phase structure in the sum over chiral blocks. 



1 Introduction 

One of the most exciting developments in the past few years has been the conjecture 
of Ooguri, Strominger, and Vafa [2, relating a suitably defined partition function of 
supersymmetric black holes to computations in topological strings. In particular, they 
argue that the perturbative corrections to the entropy of four dimensional BPS black 
holes, in a mixed ensemble, of type II theory compactified on a Calabi-Yau X are 

*E-mail: jafferis@string.harvard.edu, marsano@physics.harvard.edu 
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captured to all orders by the topological string partition function on X via a relation 
that takes the form 



Zbh ~ \Ztop\ 2 (1) 

Moreover, since the left-hand side of Q is nonperturbatively well-defined, it can 
also be viewed as providing a definition of the nonperturbative completeion of \Zt op \ 2 - 

One of the first examples of this phenomenon is the case studied in |3J |3] of type 
II " compactified" on the noncompact Calabi-Yau 0{— p) © 0{p — 2 + 2g) — * T, g , with 
S g a Riemann surface of genus g. The black holes considered there are formed from N 
D4 branes wrapping the 4-cycle O(-p) — ► S s , with chemical potentials for DO and D2 
branes turned on. The mixed entropy of this system is given by the partition function 
of topologically twisted J\f = 4 U(N) Yang Mills theory on the 4-cycle, which was 
shown to reduce to q-deformed U(N) Yang Mills on S^ 1 . This theory is defined by the 
action 

S = ^—( tr$AF+^- ( tr$AK §— / ti$ 2 AK, (2) 

9 Y m Jz a 9ym Jz a 2 9ym J^g 

where K is the Kahler form on £ 5 normalized so that E 5 has area 1, $ is defined 
to be periodic with period 2tt, N denotes the large, fixed D4-brane charge, and the 
parameters g YM ,Q are related to chemical potentials for DO-brane and D2— brane 
charges by 

<t> m = ^~ ^ 2 = 4^. (3) 
9ym 9ym 

The authors of [3] went on to verify that the partition function of this theory indeed 
admits a factorization of a form similar to Q 



Z = E E Zf^ + {t+pg s l)Zf^'(i-pg s i) (4) 
tei P,P' 

where 



xY,Q E ^ KR e' W W PR (q)W P n R (q) (5) 

zfp-{t) = {-ir\ + \ pi \zf t ^(t) 

and 

9s = 9ym t = ^-^9ym n ~ *8 (6) 

where P, P', and R are SU(oo) Young tableaux that are summed over. The extra 
sums that seem to distinguish @ from the conjectured form have been argued to be 

1 Before this connection was found, g-deformed Yang-Mills had previously been introduced, with different 
motivations, in [5]. 
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associated with noncompact moduli and are presumably absent for compact Calabi- 
Yau. The chiral blocks Zp^pf ,+ can be identified as perturbative topological string 
ampltiudes with 2 "ghost" branes inserted [3]. 

This factorization is analagous to a phenomenon that occurs in pure 2-dimensional 
U(N) Yang-Mills, whose partition function can also be written at large N as a product 
of "chiral blocks" |B] [Zj [B] that are coupled together by interactions analagous to the 
ghost brane insertions in @. Indeed, because of the apparent similarity between the 
pure and (/-deformed theories, it seems reasonable to draw upon the vast extent of 
knowledge about the former in order to gain insight into the behavior of the latter 2 . 

One particularly striking phenomenon of pure 2-dimensional Yang-Mills theory is 
a third order large N phase transition for T* g = S 2 that was first studied by Douglas 
and Kazakov pQ. The general nature of this transition is easy to understand from a 
glance at the exact solution 



In particular, we see that the system is equivalent to a " discretized" Gaussian Her- 
mitian matrix model with the rtj playing the role of "eigenvalues". Since the effective 
action for the "eigenvalues" is of order iV 2 , the partition function is well-approximated 
at large N by a minimal action configuration whose form is determined by a competi- 
tion between the attractive quadratic potential and the repulsive Vandermonde term. 
This configuration is well-known to be the Wigner semi-circle distribution and accu- 
rately captures the physics of (J7| at sufficiently small 't Hooft coupling. As the 't Hooft 
coupling is increased, however, the attractive term becomes stronger and the dominant 
distribution begins to cluster near zero with eigenvalue separation approaching the 
minimal one imposed by the "discrete" nature of the model. When these separations 
indeed become minimal, the system undergoes a transition to a phase in which the 
dominant configuration contains a fraction of eigenvalues clustered near zero at mini- 
mal separation. Sketches of the weak and strong coupling eigenvalue densities may be 
found in figure ^ 

As mentioned before, this theory admits a factorization of the form Zym = \Z+\ 2 
at infinite N and, moreover, the partition function can be reliably computed in the 
strong coupling phase as a perturbative expansion in 1/N jH] [Zl 0- In fact, roughly 
speaking the chiral blocks correspond to summing over configurations of eigenvalues to 
the right or left of the minimally spaced eigenvalues near zero. At the phase transition 
point, though, it is known that the perturbative expansion breaks down ^0] so that 
the perturbative chiral blocks cease to capture the physics of the full theory, even when 
perturbative couplings between them are included. 

A natural question to ask is whether the chiral factorization (^Q) exhibits similar 
behavior in any known examples. The goal of the present work is to lay the groundwork 
for studying this by first addressing the following fundamental questions. First, does 
a phase transition analagous to that of Douglas and Kazakov occur in the g-deformed 
theory on S 2 ? If so, is there a natural interpretation for the physics that drives it? 
Is there a correspondingly interesting phase structure of the perturbative chiral blocks 

2 A nice review of 2D Yang-Mills can be found, for instance in [S] 
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that allows one to catch a glimpse of this physics? Can we extend our study of the 
phase structure to nonzero values of the 9 angle? 

We will find that the answer to the first question is affirmative for p > 2. Moreover, 
we will demonstrate that, as in the case of pure Yang- Mills on S 2 the transition 
is triggered, from the weak coupling point of view, by instantons 3 . Using this obser- 
vation, we will study the theory at nonzero 9 angle and begin to uncover a potentially 
interesting phase structure there as well. To our knowledge, this has not yet been done 
even for pure Yang-Mills so the analysis we present for this case is also new. We then 
turn our attention to the trivial perturbative chiral block and find a phase transition 
at a value of the coupling which differs from the critical point of the full theory. In 
particular, for the case p = 3, which we analyze in greatest detail, the perturbative 
chiral block seems to pass through the transition point as the coupling is decreased 
without anything special occurring. 

The skeptical reader may wonder whether a detailed analysis is required to demon- 
strate the existence of a DK type phase transition in g-deformed Yang-Mills as it may 
seem that any sensible "deformation" of pure Yang-Mills will retain it. However, we 
point out that the g-deformation is not an easy one to "turn off". In fact, from the 
exact expression for the partition function of the g-deformed theory 



ZbH = ZqYM,g = 

(8) 



we see that, in order to make a connection with the 't Hooft limit of pure Yang- 
Mills it is necessary to take the limit gyu — > 0, N — > oo, and p — * oo, while holding 
ffy]tfNp, which becomes identified with the 't Hooft coupling of the pure theory, fixed. 
Because we are interested in finite values of p, we need to move quite far from this 
limit and whether the transition extends to this regime is indeed a nontrivial question. 
However, a glance at (JHJ gives hope for optimism as again the effective action for the 
rii also exhibits competition between an attractive quadratic potential and repulsive 
term. We will later see that this optimism is indeed well-founded. 

The outline of this paper is as follows. In section 2, we review the analysis of 
the Douglas-Kazakov phase transition in pure 2-dimensional Yang-Mills theory with 
9 = that was heuristically described above and initiate a program for studying the 
phase structure at nonzero 9. The primary purpose of this section is to establish the 
methods that will be used to analyze the g-deformed theory in a simpler and well- 
understood context though, as mentioned before, to the best of our knowledge the 
extension to nonzero 9 is novel. In section 3, we study the g-deformed theory and 
demonstrate that it indeed undergoes a phase transition for p > 2 at 9 = and make 
steps toward understanding the phase structure at nonzero 9. In section 4, we turn 
our attention to the chiral factorization of the g-deformed theory as in Q and study 
the phase structure of the trivial chiral block, namely the topological string partition 
function on O(-p) ® 0(p — 2) — > P 1 . We conjecture that, for p > 2, this quantity itself 
undergoes a phase transition but at a coupling smaller than the critical coupling of the 
full g-deformed theory. Finally, in section 5, we make some concluding remarks. 



3 The contribution of instanton sectors in pure Yang-Mills was first written in |12j 
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Weak Coupling 



n/N 



Strong Coupling 



n/N 



Figure 1: Sketch of the dominant distribution of rij in the weak and strong coupling phases 
of 2-dimensional Yang-Mills on S 2 



While this work was in progress, we learned that this subject was also under inves- 
tigation by the authors of [T2|, whose work overlaps with ours. 

2 The Douglas-Kazakov Phase Transition in Pure 
2-D Yang-Mills theory 

In this section, we review several aspects concerning two-dimensional U (N) Yang-Mills 
theory and the Douglas-Kazakov phase transition with an eye toward an analysis of the 
^-deformed theory in the next section. In addition to reviewing results for 9 = 0, we 
perform a preliminary analysis of the phase structure for nonzero 9 using an expression 
for the partition function as a sum over instanton contributions. 

2.1 Review of the Exact Solution 

We begin by reviewing the exact solution of 2-dimensional Yang-Mills. The purpose for 
this review is to obtain an expression for the partition function as a sum over instanton 
sectors that will be useful in our later analysis. There are many equivalent ways of 
formulating the theory. For us, it will be convenient to start from the action 

S=^—[ tr$AF+-?— [ tr& A K J—/ tr$ 2 AK (9) 

9ym Js 2 9ym Js 2 ^9ym Js 2 

where $ is a noncompact variable that can be integrated out to obtain the standard 
action of pure 2-dimensional Yang- Mills theory. To obtain the exact expression Q , we 
first use gauge freedom to diagonalize the matrix <£, which introduces a Fadeev-Popov 
determinant over a complex scalar's worth of modes 



A FP = det ([&*]) 



(10) 



and reduces the action to 



S'Diag * = f <t>a [dA aa - iA a/3 A Ap a ] + [ <f> a - r"4 / (H) 

9ym Jt, 9 Y m 2 9 Y m Jt, 

Integrating out the off-diagonal components of A in yields an additional de- 
terminant over a Hermitian 1-form's worth of modes 



det f F 1/2 ([<M) (12) 
which nearly cancels the determinant (|10|). Combining Q1U|) and (|12j) . we are left 



with only the zero mode contributions 

A(0) 26 ° 



A(</>) 2 (13) 



A(</>) bl 

where A(0) is the usual Vandermonde determinant 

a(^= n ( m ) 

l<i<j<N 

As a result, we obtain an Abelian theory in which each .F a is a separate i/(l) field 
strength |3| 



Z = y A(</>) 2 expjy 



S 2 , 



(15) 



To proceed beyond (|T5|) . let us focus on integration over the gauge field. The field 
strength Fi can be locally written as cL4j but this cannot be done globally unless Fi 
has trivial first Chern class. We thus organize the gauge part of the integral into a 
sum over Chern classes and integrations over connections of trivial gauge bundles. In 
particular, we write Fi = 2nriK + F[ where rj £ Z and F[ can be written as dAi for 
some Ai. In this manner, the third term of (|15|) becomes 

— — f pirnKta + fcdAi) (16) 
9 Y m 

Integrating by parts and performing the Ai integral yields a 5 function that restricts 
the (pi to constant modes on S 2 . As a result, we obtain the following expression for 
the pure Yang-Mills partition function organized as a sum over sectors with nontrivial 
field strength 

Z = Y i Z? (17) 

f 

where 

N 

Z?= / [J#iA(^) 2 exp{-^ M ^ -i&(0 + 27rrO} (18) 

i=l 
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where we have Wick rotated (for convergence) and rescaled (for convenience) the 
fa. The sum over n can now be done to yield a 5- function that sets fa = nt for integer 
Hi, leaving us with the known result for the exact 2D Yang-Mills partition function 

Z= £ AK)^)exp|-^^n?-z^nA (19) 

As noted in the introduction, this has precisely the form of a " discretized" Hermitian 
matrix model, with the fa playing the role of the eigenvalues. In particular, we note 
that the only effect of summing over in Q18|) is to discretize the eigenvalues of this 
matrix model. Because it is precisely this discreteness that will eventually give rise to 
the phase transition, it is natural to expect that the transition can also be thought of 
as being triggered by instantons and the critical point determined by studying when 
their contributions become nonnegligible. We will come back to this point later. 

Finally we note that the 9 angle has a natural interpretation as a chemical poten- 
tial for total instanton number. This is easily seen by shifting fa by %Qjg\ M in the 
contribution from a single instanton sector, to obtain 



N 

n 



' exp < 



9ym 



2tt6 



9ym 



NO 2 
%9ym 



(20) 



2.2 The Phase Transition at 6 = 

In this section, we review the derivation of the transition point in the theory with 6 = 
using an analysis that we will eventually generalize to study the g-deformed theory. 
We can write the exact result (IT9l) as 



Z= e W {-N 2 S eff ( ni )} (21) 

where 



2N 2 ^ y 3 ' 2N ^ \N 



(22) 



At large N, the sums can be computed in the saddle point approximation. Extrem- 
izing S e ff(m) yields 

ArEx-T-^. = (23) 



N fri.tjH- fa 2N 
Following standard methods, we define the density function p 



(24) 
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and approximate p(u) by a continuous function at large iV for distributed in a 
region centered on zero. In terms of p, 1(23(1 becomes 

p[dx-^ = -±u (25) 
while the fixed number of eigenvalues leads to the normalization condition 

dup{u) = l (26) 
To solve the system (j55)l . (|2^|) . it is sufficient to determine the resolvent 

v(u) = t dx (27) 

J a x ~ u 

as p(u) can be computed from v(u) by 

p(u) = — lim (v(u + ie) — v(u — ie)) (28) 

e-^o 2m 

The resolvent for the system 1)25(1 , 1)26(1 is of course well-known (see for instance )14j ) 



u(u ) = __^_y u 2_ i (29) 

and leads to the familiar Wigner semi-circle distribution 



As pointed out by Douglas and Kazakov pQ, the discreteness of the original sum 
imposes an additional constraint on the density function, p, namely 

p{u) < 1 (31) 

which is satisfied for A such that 

A < vr 2 (32) 

For A > 7T 2 . (pffll) ceases to be an acceptable saddle point of the discrete model. 
Rather, the appropriate saddle point in this regime is one which saturates the bound 
1)31(1 over a fixed interval. Douglas and Kazakov performed this analysis and found 
an expression for the strong coupling saddle in terms of elliptic integrals. Since the 
weak and strong coupling saddles agree at A = tt 2 , it is clear that the order of the phase 
transition which occurs at A = tt 2 must be at least second order. With the saddles 
in hand, it is also straightforward to compute F'(X), the derivative of the free energy 
with respect to A, as it is simply proportional to the expectation value of nf. On the 
saddle point, this becomes 

F'(A) ~ N 2 J p{u)u 2 (33) 

From this, we immediately see that continuity of the saddle point distribution 
through the critical value also implies that F'(X) is continuous and that, in fact, the 
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transition is at least of second order. To go further, it is necessary to obtain the strong 
coupling saddle and evaluate ()33j) in both phases. Douglas and Kazakov have done 
precisely this and demonstrated that the transition is actually of third order. 



2.3 Instanton Sectors and 6^0 

We now attempt to study the theory with 8^0. It is not clear to us how to im- 
pose the discreteness constraint on the saddle point distribution function p once the 
action becomes complex so we need to probe the phase structure with nonzero 8 by 
another means. The key observation that permits us to proceed is that of who 
demonstrated that the phase transition in the 8 = theory is triggered by instantons 
in the sense that it occurs precisely when instanton contributions to (|17j) are no longer 
negligible at large N. Such a result is not surprising given that it is the discreteness of 
the model, which arises from the sum over instanton sectors, that drives the transition. 

To study what happens with 8 ^ 0, we therefore turn to the instanton expansion 
(|17|) and ask at what value of the coupling does the trivial sector fail to dominate. Of 
course, we must be careful since taking 8 — > —8 is equivalent to taking rj — > — rj and 
shifting 8 —* 8+2irk is equivalent to shifting — > ri+k for integer k so, while the trivial 
sector dominates near 8 = 0, the sector (— n, —n, . . . , — n) dominates at 8 = 2irn. We 
avoid any potential ambiguities by restricting ourselves to < 8 < ir. There, we expect 
the trivial sector to dominate for A below a critical point at which the first nontrivial 
instanton sectors, corresponding to f = (±1,0, ... ,0), become nonnegligible. It is the 
curve in the X/8 plane at which this occurs that we now seek to determine. 

We begin with the trivial sector, which is simply the continuum limit of the full 
discrete model (|19j) . From ()'2()|) . the partition function within this sector is given by 



r N \ a 2 N8 2 I 

^=(o,o,..,o) = J II ex P \ E - <t>i? - 9j Y~ E # - 2^ I ( 34 ) 



i=l I i<j 



The saddle point equation for the (pi integral is dominated by the Wigner semicircle 
distribution found in the previous section (|3l)|). 

We now turn to the family of instanton sectors with r = (r, 0, . . . , 0) (we will 
eventually take r = ±1): 

N 

Zf=( r ,o,...,o) = / []^iexp{-iV 2 % (ri0r .. ]0) } (35) 

i=l 

where 



i<j i y ' K ' 

To evaluate (|35|). we note that the saddle point configuration for (fa, . . . , (f>^ in (|35|) 
is precisely the same as that for the trivial sector since the only difference between 
the effective actions in these sectors lie in O(N) out of the 0(N 2 ) terms. We may 
thus proceed by evaluating the 02, • • • > 4>N integrals using (f3"U|) and performing the 4>i 
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integral explicitly in the saddle point approximation 4 . The effective action for the (j>i 
integral becomes 

S ef f(<h.) = -J dyp(y)H ( j )1 -y) 2 + ^l-2nir(^ 1 -j^-'^- (37) 
and the saddle point value of 0i is determined by 

- 2-Kir = (38) 

For A < 7T 2 , there is one saddle point, which lies along the imaginary axis at 

2i I 

Oi)A<7r 2 r 2 = ysign(r)V r 2 vr 2 - A (39) 

As A increases, this saddle point moves toward the real axis, eventually reaching it 
at A = 7r 2 r 2 and splitting in two 

(<h)x>** =±^ V / A-vr 2 r 2 (40) 

We now wish to obtain the real part of the "free energy" in these nontrivial sectors 
by evaluating the effective action (|3*7)l on these saddle points. Comparing with the 
corresponding free energy of the trivial sector, we obtain 



N (i^=(±l,o,.,o) " ^=(0A-,0)) = — + -jp7 [y^ ) A < rV 

2ir8r , o 9 

A > 7T 7* 



(41) 



R ' 2 \l- ^T^J V 4 



A 

where j(x) is defined as in Gross and Matytsin jllj 

' , ( ' + vT^7\ n n ^ ( l-x) 2s , 

4s 2 -1 

s=l 

For 6 = 0, this result is in agreement with that of and demonstrates that the 
partition function of the one-instanton sector is exponentially damped compared to 
that of the trivial sector at large iV for A < it 2 . For A > 7r 2 , on the other hand, these 
sectors contribute with equal magnitudes at large iV and hence instantons are no longer 
negligible. 

In addition this this, we are now able to see at least part of the phase structure 
with nonzero 6 as well. While we have only considered a restricted class of instantons 
here, it is natural to assume that, as A is increased from zero, the most dominant of 
the nontrivial sectors will continue to be those of instanton number 1 and hence that 



4 Since we are only interested in the magnitude of the partition function of each sector, it will suffice 
to determine the real part of the free energy, which is obtained by simply evaluating the action at the 
potentially complex saddle point. In particular, we will not need to worry about the precise form of the 
relevant stationary phase contour. 
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they will continue to trigger the transition away from 9 = 0. We thus arrive at the 
following phase transition line in the \/9 plane, which is plotted in figure [2] 



Using the relation between shifts of 9 and shifts of the rj, we can now extend the 
phase transition line (|43|) outside of the region < 6 < tt. Moreover, we can attempt 
to guess the form of the phase diagram beyond this line. While we know that there 
are no further transitions at 6 = 0, it seems likely to us that this is not the case for 
6 7^ 0. The reason for this is that for 9^0, it is no longer true that several different 
sectors contribute equally beyond the transition as is the case for 9 = 0. Rather, 
because of the ^-dependent shift in the free energy, there will in general be one type of 
instanton sector which dominates all others for any particular value of A and, as A is 
increased, the instanton number of the dominant sector also increases. As a result it 
seems reasonable to conjecture that there are additional transitions that smooth out 
at 9 = 0. 

We can now ask when sectors of higher instanton number begin to dominate the 
partition function. If we assume that the only relevant instantons for determining the 
full phase diagram are those of the sort considered here, then it is easy to proceed. How- 
ever, we find it quite unlikely that the multi-instanton sectors that we have neglected 
in the present analysis remain unimportant throughout. To see this, consider following 
a trajectory at fixed small A > and along which 9 is increased from 0. For 9 less than 
a critical value 9 C given by (|43|) . the dominant sector is the trivial one. Just beyond 
this value, the dominant sector is that with f = (—1, 0, . . . , 0). Proceeding further, we 
expect to hit additional critical points at which a sectors with larger instanton number 
begin to dominate. Eventually, though, we will approach 9 = 2ir — 9 C , beyond which 
the (—1, —1, . . . , —1) sector becomes the dominant one. More generally, for each region 
in the range < 9 < it in which an instanton sector f dominates, there is a region in 
the range tt < 9 < 2ir in which an instanton sector f = (l,l,...,l) — r dominates. A 
conjecture consistent with this is that a series of transitions occur along this trajectory 
in which the dominant sectors take the form (—1, —1, . . . , —1, 0, 0, ... , 0). We sketch 
a phase diagram based on this conjecture in figure |21 A deeper analysis which takes 
into account the multi-instanton sectors not considered here in order to confirm/reject 
a picture of this sort would be very interesting, but is beyond the scope of the present 
paper. 



We now proceed to study the g-deformed theory gous to the previous 

section 5 . Once again, we begin by reviewing the exact solution while obtaining an 
expression for the partition function as a sum over instanton sectors. 

5 In particular, we will once again use saddle points to study the large N behavior of various instanton 
sectors. In pure Yang-Mills, one can obtain more precise results by the technique of orthogonal polynomials 
This can in principle be generalized to the g-deformed case due to the recent identification of an 
appropriate set of orthogonal polynomials for such an analaysis |15| 




(43) 



3 The deformed Theory 



n 




(-l,0,-,0) Dominates 
1 



(-1 ,-l,0,-,0) Dominates? 

2 



(0,0,. ..,0) Dominates 



(-1,..., -1,0,0) Dominates? 

(-1,...,-1,0) Dominates 




(-1,-1,...,-!) Dominates 



Figure 3: Conjecture for the 2D Yang-Mills phase diagram provided we take the curve 
seriously over the entire range < 9 < tt 2 . The dotted lines represent the first of numerous 
conjectured transition lines that cluster near 9 = n. 
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3.1 Review of the Exact Solution 

We begin by recalling the action of the g-deformed theory 
1 /"...„ 9 /"...„ p 



S 



tr $ A F + 



tr $ A if 



tr$ 2 A if 



(44) 



ffy-M Js 2 9ym Js 2 2g YM Js 2 

where <3? is periodic with period 2tt. As demonstrated in the analysis proceeds 
exactly as in section 2.1 with the periodicity of giving rise to the replacement 



|e i(<fc-^)/2 _ e -i(^-^)/2 



a(^) - A(^) = n [. 

l<i<j<N 

In particular, (|17 |) -l|18 |) become, after an identical Wick rotation and rescaling 



(45) 



(46) 



with 



JqYM,r 



N 

n 

i=l 



n > 

l<i<j<N 



-QYM^-'t'j)/ 2 — e -9YM < y't>3- ( t > i)/ 2 



x™P< l - i ^<fi-i(e + 2irr i )<j> i } (47) 



Performing the sum over the restricts 4>i to integer values, as in the case of pure 
Yang-Mills, and yields the known result 



Z q YM = ^ 



l<i<j<N 



xexpj-^J> 2 -z0j>l (48) 



As before the partition function resembles a matrix model, this time what appears 
to be a "discretized" Hermitian matrix model with unitary measure 6 . In addition, 
the sum in Q49|) can again be interpreted as a sum over instanton sectors with the 
trivial sector corresponding to the continuum limit and capturing all aspects of the full 
partition function except for the discreteness. Moreover, by shifting <pi by i9/g YM p 
in (|49|) . we see that the 9 angle continues to carry the interpretation of a chemical 



6 Matrix models of this sort have been studied before in the context of topological strings, first in J|j] and 
later in J7]. Their "discretization" has also previously been studied in 
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potential for total instanton number: 



J qYM,r 



N 

n 



l<i<j<N 



x exp < 



9ymP 



2tt6 



2 5y M P J 



(49) 



It seems reasonable to expect that the trivial sector will be dominated by a saddle 
point configuration analagous to the Wigner semicircle distribution for pure Yang-Mills 
below a critical coupling. Beyond this point, one might expect the attractive term in 
the potential to become sufficiently large that this distribution is too highly peaked 
to be consistent with discreteness, at which point we expect to find a phase transition 
that, due to its connection with discreteness, can again be thought of as arising from 
the effects of instantons. 



3.2 The Phase Transition at 6 = 



We now specialize to the case 9 = for simplicity and determine the saddle point 
configuration for small values of the coupling. We will find a distribution analagous 
to the Wigner semicircle distribution which will, for sufficiently large coupling, violate 
the constraint arising from discreteness leading to a phase transition analagous to that 
of Douglas and Kazakov. To proceed, we write the exact result (|48jl as 



Z= £ exp{-N 2 S eff (n t )} 



(50) 



where 



S, 



eff 



2N 2 



In ( 2 sinh 



\{rii 



2N 



+ O AT 2-^1 \ N 



(51) 



We proceed to study this in the saddle point approximation. The saddle point 
condition is easily obtained 



1 

N 



coth 



X(m - rij) 



2N 



N 



(52) 



We again follow standard techniques and introduce a density function p(u) accord- 
ing to (|24|) in terms of which (|52|) can be written 



P f dy p(y) coth bt_A 
and which must satisfy the normalization condition 

dyp{y) = 1 



px 



I 



(53) 



(54) 
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To solve (|53[) and (|54|). it is convenient to make the change of variables Y = e Xy . 
In this manner, we may rewrite the system (|53|) . (|54|1 as 



P / dY 



dY 



Y -X 

m 



2A 



In 



Xe~ x ' p 



(55) 



where 



~p{Y)=~p(e x y) = \p(y) 



(56) 



As usual, solving this system of integral equations is equivalent to determining the 
resolvent 



v(X) 



dY 



P(Y) 



(57) 



Y-X 

Fortunately, this model has been studied before and the resolvent found to be 



v(X) = y In 
A 



X + l + ^jX -a+)(X -q~ 
2X 



where 



a± = -l + 2e x/p l±Vi 



-X/p 



(58) 



(59) 



From this, we find that p(u) is nonzero only for a_ < e Xy < a+, where it takes the 
value 



p(y) = 7T~ arccos 

Z7T 



(60) 



The appropriate branch of the arccos function for this solution is that for which 
< arccos x < tt. As in the case of pure Yang-Mills, the discreteness of the model 
imposes an additional constraint 



p{u) < 1 (61) 

which is satisfied provided 

A < X cr u = —pin cos 2 — (62) 
P 

For A > Xcrit, (j6*U)) ceases to be an acceptable saddle point of the discrete model 
and we expect a phase transition to a distribution analagous to the strong coupling 
saddle of Douglas and Kazakov at this point. Note that X cr it is finite and nonzero only 
for p > 2 so we conclude that a phase transition of the Douglas-Kazakov type occurs 
only for these p in the g-deformed model. 

To go beyond X cr u, we must look for a new saddle which saturates the bound (l6*Tj) 
over a fixed interval. As in the pure Yang-Mills case, this strong coupling saddle must 
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agree with that at weak coupling at A = A cr ^. Moreover, we note that the derivative 
of the free energy with respect to A is given at large TV by 



F'(X) = N 2 

N' l p 



1 f f \(x — y)\ Xp f 
-P / dxdy p(x)p(y)(x - y) coth ( - )+y / d V pWlV 1 

dyp{y)y 2 

(63) 



where we have used ()53(l . Consequently, we see that, as in pure Yang-Mills theory, 
continuity of p(y) through the transition point guarantees that F'(X) is continuous 
and hence the transition must be of at least second order. It seems very likely to us, 
however, that the transition will continue to be of third order as in the case of pure 
Yang-Mills. We include a brief discussion of the necessary tools to study the strong 
coupling saddle point in Appendix A. A more complete presentation, as well as further 
arguments for the third order nature of the transition, can be found in jl.Sj . 



3.3 Instanton Sectors and 6 ^ 

We now turn 9 back on and attempt to study the phase structure. As in the pure 
Yang-Mills case, we seek to probe the phase structure by analyzing the instanton 
expansion (f4*9|) . Again, the trivial sector dominates at small 9, X and we expect the 
phase transition to occur when the first nontrivial instanton sector, corresponding to 
r = (±1, 0, . . . , 0), makes a nonnegligible contribution to the partition function at large 
N. We now work by analogy and seek to determine the curve in the X/9 plane that 
bounds the region in which the trivial sector dominates. 

We begin with the trivial sector, which is simply the continuum limit of the full 
discrete theory (|4*5)l 



V=(0,0,...,0) 
N 



[ Y[ d(f>i exp I ]- ^2 ln ( 2 sinh 
3 i=i \ i*i V 



X((f)i - <fij 
2N 



Xp 
2N 



E 



N 



N9 2 \ 
2 9ymP J 



(64) 



The saddle point distribution for this integral is given by the function p found in 
the previous section (|6U|) . 

We now turn to the family of nontrivial instanton sectors with r = (r, 0, . . . , 0) (we 
will later take r = ±1) whose partition functions take the form 



J r=(r,0,...,0) 



N 

/ 11 dxj>i exp{-N 2 S, 
J 1=1 



e//,M,...,0) 



} 



(65) 



with 



S, 



e//,(r,0,...,0) 



— ^ In sinh 2 
N 2 ^ 



i<j 



X((pi - <f>j) 



N 



Xp 2 27rir 



i=l 



Xp 



2Xp 
(66) 
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The integrals over </>j with i > 1 are dominated by the saddle point (|6Uj) leading to 
the following effective action for the integral 



?eff{9h 



J dyp(y) 



In sinh 2 



A(0i - y) 



+ -J-</>i - 27r«r 



i0\ 



A^ 2 



2 J ' 2 T1 "" V Ti Apy 2Ap 
It is easy to determine the saddle point equation for <f>\ using (|53|) 



A 
P 



1 - 2v ( e^ 1 



2ixir 
P 



(67) 



(68) 



where t> is the resolvent from (|57j). Solving this equation is straightforward, though 
requires a bit of care. If we write r = np + m where n, m are integers and m £ 
(— p/2,p/2], then the solution can be written as 



2» . / \ 
?i = — sign(r) arctan 
A 



-A/p 



scc- 



irm 



P 



+ 



2mn 



A 



1±W1 



-A/p 



7rm 



sec z 



+ 



e - A /Psec 2 
2-Trin 



Trm 



> 1 



-A/p 



7T77T, 



A 



sec 



P 



< 1 
(69) 

where the appropriate branch of arctan in (|69|) is —tt/2 < arctan x < ir/2 and that 
of In is such that it has zero imaginary part for real argument. 

As in the case of pure Yang-Mills, the saddle points lies along the imaginary axis for 
small A and, as A approaches the critical point \ cr it = —pin cos 2 ^, the saddle points 
corresponding to r = ±1 hit the origin and splits in two along the real axis as A is 
increased further. 

To study when particular instanton sectors become nonnegligible, we must evaluate 
the (real part of) the free energy in each sector, namely 



F*= Re 



2iTr9 



where 



dy p(y) In sinh 2 



2 / X(x - y) 



(70) 



(71) 



In Appendix B, we obtain an exact expression for the function g{4>) and, in addition, 
simple integral formulae for the real part of g((p) for <f> along the real and imaginary 
axes. Using the results of Appendix B, we may write the difference between the real 
parts of free energies of the r = (r, 0, . . . , 0) and r = (0, 0, . . . , 0) sectors as 



AT 



f=(r,Q,...,0) 



F, 



r=(0,0, 



.,0)) 



4p 



on 



da 



irr 



sign(r) arccos ( e 



-\/2p\ 



for e A ^ p cos 



cos a\ 
2 irm 



+ —— (2-Knp - 9) 
Xp 



< 1 



2irr 

— — (2Trnp 
Xp 



P 

irm 



for e VP cos 2^ >1 
P 



(72) 
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where 



«o = sign(r) arccos ( e 



,\/2p 



ixm 

cos 

V 



(73) 



and the appropriate branch of arccos is such that < arccos x < ir. 

Let us first look at the case 9 = 0. The integrand in 1)72(1 is manifestly positive 
definite for ao nonzero and e A//p cos 2 (7rm/p) < 1 so the integral is positive definite for 
ao > and for ao = 0. As a result, we see that as the coupling is increased from 
zero, the trivial sector dominates until the critical point A = X cr u ITS is reached, at 
which point ao becomes zero for r = ±1 and the difference in free energies vanishes for 
the one-instanton sectors r = (±1, 0, . . . , 0). Consequently we conclude that, as in the 
case of pure Yang- Mills, the phase transition at A = A cr it (|62|) is triggered by nontrivial 
instanton sectors becoming nonnegligible and contributing with equal magnitude to 
the trivial one at large N. 

Now, what happens when we turn 9 back on? It is again reasonable to assume 
that, as A is increased from zero, multi-instanton sectors of the type not considered 
here remain unimportant so that most dominant of the nontrivial sectors for A < A cr it 
(|62(1 will continue to be those of instanton number ± 1. This leads us to the following 
phase transition line in the X/9 plane, which is plotted for p = 3 in figure |I] 7 



2p 2 f a ° 

da 



TT TT 2 Jo 



sign(r) arccos I e " y \ cos a 

p ' 



-V2PI 



(74) 



The qualitative behavior of this curve is identical to that of the pure Yang-Mills 
case in that the critical coupling decreases from X cr u Ijfi2|) at 9 = to zero at 9 = 
tt. Following the same reasoning applied to pure Yang-Mills, we also expect further 
transitions beyond the line (|74j) beyond which sectors with larger instanton number 
dominate the full partition function. The required symmetries under simultanoues 
shifts/flips of 9 and r then lead us to conjecture a phase diagram of the form © in 
the case of the g-deformed theory as well. It would be interesting to pursue a further 
analysis in order to verify or reject this picture. 



4 Chiral decomposition, and phase transition in 
topological string 

One of the most interesting features of the (/-deformed Yang- Mills theory we have been 
examining is its factorization into chiral blocks at large N. This occurs in an analogous 
manner to the chiral decomposition of pure Yang Mills, with the new ingredient being 
its identification by 3 with the topological string partition function. We will show 
that the chiral partition function has a natural interpretation in the strong coupling 
phase as a limiting description of the eigenvalues on a single side of the clumped region. 

In contrast to the case of Yang-Mills on a torus, where a similar picture was devel- 
oped [1] ^S], here of course the Fermi surfaces do not decouple, even perturbatively, 
because of the nontrivial repulsion between the eigenvalues. We will see that this is 
encoded in the sum over " ghost" branes found in [3] . The presence of the weak coupling 

We would like to thank M. Marino for correspondence regarding a mistake in the plotted curve in version 
1, which was the result of a sign error. 
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Figure 4: Plot of the curve for r = l,p = 3. 



phase adds a new feature to this picture, since at some value of the 't Hooft coupling, 
the two Fermi surfaces simply vanish. We will investigate what effect this has on the 
trivial chiral block. Throughout this entire section, we consider only the case 8 = for 
simplicity. 

We now look a bit closer at the trivial chiral block, Z qYM,+ , which is to be identified 
with the topological string partition function without any ghost brane insertions 3 . 
Using ®, we see that it takes a particularly simple form 



R 



Y^q^e-^Wniq- 1 )' 



(75) 



R 



where the sum is over Young diagrams, R and the objects appearing in the sum are 
defined as follows: 



KR = Y J R i( R i~ 2i + 1 ) 

i 

W R {q- 1 ) = s R (q*-^ 



(76) 



q = e 



9ym 



where Ri denotes the number of boxes in the ith row of the diagram R and s^x 1 ) 
is the SU(oo) Schur function associated to the diagram R. 

It is clear that, if we include all Young diagrams, the series (f?3|) diverges for any q 
for p > 3, since arbitrarily negative powers of q appear for diagrams with sufficiently 
large diagrams. This should not disturb us, however, because the topological string 
is merely an asymptotic expansion, and only the perturbation theory about g s = is 
well defined. Thus we truncate the sum over Young diagrams to chiral representations 
of U(M), for some finite M, and take a 't Hooft limit keeping A = g s M fixed. In fact, 



19 



the chiral decomposition (jHJ) indicates precisely what this cutoff should be as one can 
easily verify that @ holds at finite N only when the Young diagrams R are taken to 
have at most N/2 rows. To be sure, the model that we study is the chiral block defined 
in this way. We then assume that, based on (J3J), this is the proper manner in which to 
define the perturbative topological string partition function. 



4.1 Interpretation of the trivial chiral block 

Heuristically, we expect the chiral partition function to describe a system in which only 
the eigenvalues on one side of the clump in the strong coupling phase are " dynamical" 
in the sense their locations are summed. We can isolate such contributions by placing 
an infinite clump of eigenvalues about the origin which drives the two fermi surfaces 
apart, effectively decoupling them. Of course, one must be careful when doing this to 
avoid introducing divergences in the partition function. A particularly natural way to 
proceed is to start with the (/-deformed Yang-Mills partition function expressed as a 
sum over U(N) Young diagrams. This is easily obtained from the exact expression (JHJ 
by the change of variables 



TZi = rii- pi (77) 

with p the Weyl vector of U(N) 

p i = l(N-2i + l) (78) 

Note that the TZi parametrize deformations from the maximally clumped rectan- 
gular distribution so organizing the partition function as a sum over the TZi is quite 
natural for the problem at hand. The resulting expression is 

2 q ic 2 {n) e -ie\n\ (79) 

where the sum is now over U(N) Young diagrams 7Z. The quadratic Casimir and 
1Z\ are defined as 



ZqYM ~ ^2 n 



[Hi - Kj + j 
[j - 4 



C 2 (Tl) =K n + N\K\ 

\K\ =J> ( 8 °) 

i 

and [x] q is the q- analogue defined by 

[x] q = q x ' 2 - q~ x / 2 (81) 
To add an additional M eigenvalues, we simply pass from N to M + iV in the above 
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expression, obtaining a sum over U (M + N) representations that becomes 



n 

= £ 



n 

l<i<i<M+JV 



- 72.J + j - i] 



[j ~ i\ 



MM+N]\n+K U ) e -ie\n\ 



-(N+M-1)\K\ | J 

l<i<j<M+iV 



* rP^"j 3 



1-3 



Z([M+N]\K\+K Re -ie\K\ (82) 



,(M+A0 
'7J 



2 ^(M+jV)|7e|+f« We -i6W 



If we now restrict to representations with at most N/2 rows, corresponding to 
treating only the largest N/2 eigenvalues as "dynamical", and take the limit M — > oo, 
which is well-defined only for p > 2, we obtain the result 



pK R /2 



(83) 



which is our old friend (|75|) . We thus see the sense in which the chiral block describes 
the dynamics of eigenvalues to one side of the clump. 



4.2 Phase Structure of the Trivial Chiral Block 

We now perform a preliminary study of the phase structure of the trivial chiral block. 
As mentioned before, this quantity is defined with a cutoff on the sum over represen- 
tations to those with at most N/2 rows so, to study this model as g s — * 0, we must 
work in a 't Hooft limit where gy^N ~ t/2 is held fixed while taking N ~ t/2g s to oo. 
It is convenient for this analysis to return to the ni variables, of which there are now 
only N/2. The effective action in these variables becomes 



C 1 N ' 2 9 

N 2 2N 2 2^ m \ e ) + 

N/2 N/2 N/2 

^£M^) + i£Gr+!£! («> 

i=l i=l i=l 

For simplicity, we now set 9 = 0. Proceeding via the usual saddle point method, 
we find that the configuration which extremizes the action satisfies 

i=i 

Introducing the eigenvalue distribution p as usual, we have 

dyp(y) + i ln fi_ e -A,U^ = 



1 _ e X(x-y) \ \ J 2 
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where, because the cutoff is at N/2 rather than N, the normalization condition 
becomes 

J dyp{y)= l - (87) 
To solve this, we make the change of variables 



U = e Xu 



If we then define 



p(U) = jp(u = jlnuY (89) 
then the saddle point equation that we must solve becomes 

with the constraint 

» = 1 (91) 
Because we expect clumping near zero, we look for a saddle point of the form 

1 < u < a 

p(u) = Tp{u) a<u<b (92) 
u > b 

We now define the function ip as 

$(U) = ( u = jlnUj (93) 
and note that ()9U|) and (|91() may be written in terms of ifi as 

P r^n__^ u+lHS _ A) 



and 

A 



dV^iV) 2, ( 



where 

A = e Xa B = e Xb (96) 
As usual, we may determine ip(U) by first computing the resolvent 

f(U)=[dVp^, (97) 
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The expression for f(U) is given by the well-known "one-cut" solution (see, for 
instance, [Tl]): 



B) 



dS- 



(98) 



'[AM (U-S)y/(S-A)(S-B) 

with integration contour going counterclockwise around the cut [A, B]. To evaluate 
this, we move the contour away from the cut, and pick up instead the branch cuts of 
the logarithms along [—00, 0] and [—00, A], as well as the pole at S = U. The resulting 
integrals can be easily evaluated, and we find that 



W) = -ln 



U-A 



A 



P 



In 



U J A A 
/ y /(U-A) + ^(Tr- 




^B{U — A) + \J A(XJ — B) 

y/B{U-A)-y/A{U-B) 



x-y y /( U -A)-yW^B) 

Using this, we may determine ip(U) and, from that, the density p{u): 



(99) 



p(u) 



1 







p arctan 



B-e Xu 
e Xu -A 



(p — 2) arctan 



A(B—e Xu ) 
B{e Xu -A) 



< e Xu < A \ 

A<e Xu <B 

e Xu > B ) 

(100) 

where we have returned to the original variables. 

It remains to fix the position of the cut, A and B. This is done by noting that 
the normalization constraint implies /(0) = —1 + fin A while the definition of f(U) 
implies that f(U) ~ 0(U~ 1 ) as U — ► 00. Let's look at the first condition, namely 
f{U) = 0. In the limit U -»■ 0, f(U) becomes 



On the other hand, in the limit U — > 00, f(U) becomes 



f(U - 00) = P - In ( ^— - ) + ^— 1 In I v " ' v " I (102) 



B-A\ p- 2, /v^B + v 7 ^' 



A V 4 / A ^v^-v 7 ^ 
Defining 

^/B-^fA *jB + ^/A , . 

* = 2 y = - 2 , (103) 

for convenience, we obtain the polynomial equations 

x 2( P -2) (y2 _ x 2 ) = e -A/2 (1Q4) 

xy p - x = 1 (105) 
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p 



These can be solved easily in principle. As an example, let us consider the case 
= 3. From (|1U5|) we find Y = \j\fX. We now plug into ()1U4[) and look for real 



solutions such that A 
real, namely 



(X- 1 ' 2 - X) 2 > 1. There are four solutions, two of which are 



X 



v 7 ^ 



1 ± 



3/2 



where 



4 (2 



a \6 



1/3 



e- x ' 2 + 



a 



2V332/3 



a = \ 9 + v / 3e" 3A/2 V / 27e 3A - 256e 3A /2 
These solutions are real provided 

'256 



1/3 



A > (2/3) In 
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1.499 



(106) 

(107) 
(108) 

(109) 



but in this range the plus sign leads to A < 1 so it is only the minus solution that 
will interest us. Moreover, we find that at A = 1 at 



A as 1.74174 



;no) 



and hence the density (|100|) ceases to be an acceptable saddle point of the model 
below this value of the coupling. This leads us to conjecture that the p = 3 chiral 
block, and presumably the chiral blocks for larger values of p as well, undergoes a 
phase transition at a critical value of A. Note that, interestingly, this occurs far below 
the phase transition of the full (/-deformed theory, which for p = 3 corresponds to 
Acrit = 3 In 4 ~ 4.159. We plot p(u) in this case for various values of A in figure El 
Further study of the phase structure here would be very interesting. We postpone this 
to future work. 



5 Concluding Remarks 

We have found a phase transition in g-deformed Yang- Mills theory on S 2 in the 't Hooft 
limit, driven by instantons. Using saddle point techniques we evaluated the free energy 
in the weak coupling phase. We studied the structure of the contributions of various 
instanton sectors at strong coupling, and find evidence of a rich phase structure when 
a 6 angle is turned on. These phase transitions have very interesting consequences for 
the chiral factorization into a sum over topological string expansions. We discover a 
phase transition in the trivial chiral block at a different value of the coupling, implying 
a nontrivial structure in the breakdown of the factorization as one moves into the weak 
coupling phase. 

While in this work we have focused on p > 2, our analysis also clarifies certain 
aspects of the OSV relation for p = 1 and 2 as well. In those cases, the q-deformed 
Yang-Mills remains in the weak coupling phase for all A, thus the chiral decomposition 
appears suscipious. In particular, we identify the absence of a clumped region of 



24 



Figure 5: The saddle point density p(u) as a function of A in the case p = 3 for A >^ 1.74174. 



eigenvalues with the fact that the attractor value of the Kahler class of the S 2 is 
negative for p = 1. In [2], it is argued that one can perform a flop, which essentially 
turns on a 9 angle of —it, in order to make sense of this case. This is yet another 
motivation to understand the phase structure at nonzero 9 angle better. On the other 
hand, for p = 2, the attractor size is which is consistent with the fact that only in 
the limit A — > oo do we saturate the discreteness condition. 

There are several intriguing directions which deserve further research. It would 
be very interesting to fully explore the phase structure of the g-deformed Yang-Mills 
for non-zero 6, extending our analysis of the instanton sectors deeper into the strong 
coupling regime where multi-instanton corrections may become important. This seems 
particularly promising since we have seen some indictations of a rich phase diagram. 
The manner in which the chiral decomposition breaks down in the weak coupling phase 
should also be further elucidated, since we found a phase transition in the topologi- 
cal string partition functions in the trivial sector only after the coupling had moved 
well inside the regime where we expect the decomposition into chiral blocks to break 
down. We believe that this phenomenon will be found to be driven by putatively 
nonperturbative corrections to the factorization becoming important. 
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Appendix A. Two cut solution in the strong cou- 
pling phase 

We have seen that for A > \ cr it, the instanton sectors are no longer suppressed, and 
we have a transition to a phase in which some of the eigenvalues are clumped together 
at the maximum density determined by the constraint (|31|) . Following Douglas and 
Kazakov pQ, we will substitute the ansatz, 

p(u) = 1 for — a < u < a, and p(u) otherwise, (HI) 

for the eigenvalue density in the strong coupling regime. 

Plugging this into the effective action for q-deformed Yang Mills (|51|) . we find the 
new effective action 



«)) 



S eff = Y / dup{u) u 2 + y d V V 2 ~ 2 i dup(u) / dvp{v) log ( sinh 2 (-(u 

(H2) 

dup(u) J dy log ^sinh 2 (^(u - y))^j dx J d V lo § ^sinh 2 (^(x - y))^ , 

where we have set 6 = at present for convenience. The saddle point equation is 
therefore 

W' new (u) = pXu - log = \ V J dy p{y) COth Q (u " y) ) ■ (U3) 

This is a system of repulsive eigenvalues in the effective potential, W new , which has 
two wells, one on each side of the maximally clumped region, that the eigenvalues can 
fall into. Thus we will look for a symmetric two cut solution. It will be easier to first 
change variables to U = e Xu and A = e Xa . This implies that 

plog[/-log -r-i r-T = -W — p(Y)— 



1 1 1111 o' / 1/" r\," /l 1 115 

J75^-2 - u~*A* J ) 2 J Y U2Y-2-U-2Y2 
which can be simplified to 

plogU _ 210g (^) + A( i_ 2o) . P /^. m 

Defining the resolvent by v(U) = J p< ^}^ , as in equation (|57|). we are left with 
a standard problem of two cut matrix model, see for instance |14j . The solution is 
determined by solving the Riemann-Hilbert problem for the resolvent, which is char- 
acterized by falling off as U^ 1 at infinity, having period 2ixi about the two cuts, and 
taking the form v(U + ie) — v(U — ie) = 2mp(U) along the cuts. Using well known 
techniques in complex analysis, one finds that 

v{U) = 

dz plogz - 2 log (^f) + Kh ~ 2a ) I {A - U){A~^ - U)(B - U^B- 1 - U) 
2^1 U-z y (A-z)(A- l -z)(B-z)(B- 1 -z) ' 

(115) 
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where the contour encircles the cuts from (-B -1 , A^ 1 ) and (A,B). Moving the contour 
toward infinity, we pick up instead the branch cut of the logarithm along (A -1 , A), the 
pole at z = U, and the branch cut of the log U along the negative real axis. 
Putting all of this together, we obtain 

U-A^X 1 f A y/(A - U^A- 1 - U)(B - U)(B^ 1 - U) 



v(U) = plog 17-2 log ( ^zV) +A(^-2a)+2 ds 

f-0 

-py/(A - U^A- 1 -U)(B- U)(B- ! - U) / 

7— c 



{U - s)y/(A - s)(A~ 1 - s)(B - s^B- 1 - s) 



ds- 



(U - - s)(A~i - s)(B - s)(B~i - s) ' 

(116) 

This last two integrals can be expressed as elliptic integrals of the third kind. The size of 
the cut can be determined by requiring that we reproduce the correct limit as U — > oo, 
where we must have v(U) ~ and as U — > 0, where v(U) = J ^y-p(Y) = A — 2a. 

Hence we can now compute the free energy in the strong coupling phase. 
Define the elliptic integral, /, as follows 



r(>) _ , dt l(A-U)(A-i-U)(B-U)(B-i-U) 



i A -i U-t\ (A- t)(A- 1 - t)(B - t)(B~ l - t) 

2% 



where 





-U)(l 


- BU) 




-U)(l 


-AU) 


/<*- 


-A)(s- 





(117) 



{(1- AU)F[z\m] + (A 2 - l)U[n, z\m]} 



, , { AB-\ \ 2 (B-A-^A-U) 
z = sin 4 / -— ' , m = — — , and n - 



[B-A-^iA-sY \B-AJ ' (B - A)(A' 1 -U)' 

Imposing the limits derived above implies that 



f ds _ P r ds 

Ja-i yJ{A - s){A^ - s){B - s ){B^^s) ~ 2 Loo ^JlT^^^^^^^^' 

(118) 

and 



7— c 



ds 

Sy /(A - s^A- 1 ^sj(B - s)(B- 1 - s) 

sy/(A - s)(A~ 1 -s)(B- iJpR - s) ' 

which can be solved for A and B in terms of A. Note that the divergence of the 
integral up to — e exactly cancels the loge. 



2 r . (119) 

7a-i sJ(A - s)(A- 1 -s)(B- s)(B- 1 - s) 



Appendix B. The Function g(ot) 

In this appendix, we study the function 

g(a) = J dyp(y)lnsmh 2 ( ^ {a ~ y) ^ (120) 
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where p(y) is the density function ((50)) in the strong coupling phase of the q- 
deformed Yang-Mills theory at 9 = 0. To start, we note that 



g'{a) = X dyp{y) coth 



A(a - y) 



121) 



A 1 - 2v e 



where v is the resolvent (|57|) first obtained in ^7j- If we define ^4 = e Xa , then we 
can express g(a) as a function of A. Moreover, we have that 



* (A) -I (1-2*1)) 
which can be integrated in order to obtain g(^4): 



(122) 



9 A 



where 



In ( 1 - e" A/p 



In ,4 + [lnn(,4)] 2 - ^ [lnu(A)} 
A A 

^Li 2 (1 - u(A)) - ^-Li 2 (e~ x K(. \ ) ) + C 



u(A)) 



(123) 



u(A) 



1 + A + y/(A + l) 2 - 4e x /PA 



2A 



(124) 



and C is an a-independent constant. The result (|123|) is not very useful for learning 
about free energies. However, it is useful for noting that 



2m\ . . 27tz 
9 [ a + — J = 5(a) + — 



(125) 



and hence that shifting a in this manner has no effect on the real part of g(a). This 
fact is useful for the analysis described in the main text. 

What we really need in order to pursue the analysis within the text is not necessarily 
g(a) itself, but rather simply g(a) evaluated for a purely real or purely imaginary. 
Writing g as g(z = x + iy) = u(x, y) + iv(x, y) we note that 



g'{z) = d z g(z = x + iy) = d x u(x, y) - id y u(x, y) 



and hence that 



u(x,0) = J Reg\x,y = 0)+D u(Q,y) = -f lmg'(x = 0, y) + D' 



(126) 



(127) 



where D, D' are constants independent of x and y that depend only on A and p. 
Computing u(x, 0) is particularly easy 



u(x, 0) 



Xpx 2 



(128) 
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v \e" Xa ) -v e~ iXa 



For u(0, y), we have the integral 

u(0, a) = — A / tin 
which simplifies to 



Xpa 2 , 
u{0, a) = h 2p J da arctan 



e A/p gec 2 

2 



(129) 



(130) 
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